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Abstract 

Considering one of the fundamental notions of Bers' theory of pseu- 
doanalytic functions the generating pair via an intertwining relation we 
introduce its generalization for biquaternionic equations corresponding to 
different first-order systems of mathematical physics with variable coef- 
ficients. We show that the knowledge of a generating set of solutions of 
a system allows one to obtain its different form analogous to the com- 
plex equation describing pseudoanalytic functions of the second kind and 
opens the way for new results and applications of pseudoanalytic function 
theory. As one of the examples the Maxwell system for an inhomogeneous 
medium is considered, and as one of the consequences of the introduced 
approach we find a relation between the time-dependent one-dimensional 
Maxwell system and hyperbolic pseudoanalytic functions and obtain an 
infinite system of solutions of the Maxwell system. Other considered ex- 
amples are the system describing force-free magnetic fields and the Dirac 
system from relativistic quantum mechanics. 

1 Introduction 

Bers' theory of pseudoanalytic functions mainly created in fifties of the last 
century [2] offers interesting and still not fully explored tools for studying and 
solving linear elliptic equations in the plane. Recent advances in the theory and 
its applications (llj show that some abstract constructions proposed by Bers can 
be made completely explicit and applicable to important equations of mathe- 
matical physics. For example, pseudoanalytic formal powers introduced and 
studied by Bers and later on by many other mathematicians had been obtained 



only in some very special situations which represented a considerable obstacle 
for a further development of pseudoanalytic function theory. This obstacle has 
been substantially diminished in a recent work reported in [llj due to the fact 
that there has been found a method for constructing formal powers explicitly in 
a much more general situation. Moreover, it was shown that this construction 
offers a tool for calculating explicitly complete systems of solutions of linear 
elliptic second-order equations in the plane. This progress together with some 
other recent developments posed more open questions concerning pseudoanalytic 
function theory, its generalizations and applications. For example, a possibil- 
ity to develop a hyperbolic pseudoanalytic function theory with applications to 
hyperbolic equations of mathematical physics was explored in [T2] (see also [TT] 
and [S]). In the present work we use some of the results of [T2] for obtaining 
an infinite system of solutions of the one-dimensional time-dependent Maxwell 
system. Another related important open question is the development of pseu- 
doanalytic function theory in higher dimensions. This was analized in a number 
of papers (see [2], [TU], [I] and [H]), and it is clear that in this direction the 
development is barely starting. 

The aim of the present paper is to show that some fundamental ideas of 
pseudoanalytic function theory are valid in a general situation and applicable 
to linear systems of mathematical physics both elliptic and hyperbolic. The 
starting point of Bers' theory is the concept of a generating pair which in a 
sense means a substitution of a pair of "rectilinear elements" of the plane 1 and 
i by a pair of quite arbitrary "curvilinear elements" - a pair of complex functions 
F and G which only should enjoy the property of independence in the sense that 
any complex function w can be represented in the form w = ipF + tpG where 
(p and are real-valued functions. Beginning with generalizations of the first 
definitions from analytic function theory like the derivative, Bers shows that 
behind a generating pair there is always a corresponding generalized Cauchy- 
Riemann system which is usually called the Vekua or Carleman-Vekua equation. 
The knowledge of a generating pair for a Vekua equation allows one to represent 
it in another form which is the equation for pseudoanalytic functions of the 
second kind. This form is very convenient for introducing a simple formula for 
calculation of the [F, G')-derivative in the sense of Bers and of the corresponding 
antiderivative and in fact represents a cornerstone of all further constructions 
of pseudoanalytic function theory including formal powers. 

In this work we show that the concept of a generating set of functions in 
the sense of Bers is much more universal and can be introduced in relation 
with first-order systems of mathematical physics with the aid of hypercomplex 
algebraic tools. This allows one to obtain another form of a corresponding 
system analogous to that describing pseudoanalytic functions of the second kind. 
In order to generalize the concept of a generating pair it resulted to be fruitful 
to develop a slightly different approach to this concept via a certain intertwining 
relation. We consider it in the next section. Another tool implemented in this 
paper is the algebra of biquaternions. The related notations are introduced also 
in section [21 

We consider only few examples in this paper chosen in such a way that from 
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one side it becomes clear that our approach is general and is applicable to a wide 
variety of systems and from the other it can be seen that for each particular 
physically meaningful system some special interesting phenomena may occur. 
Thus, we consider the time-dependent Maxwell system for inhomogeneous media 
(section 3), the system describing force-free magnetic fields (section 4) and the 
Dirac equation (section 5). In the case of the Maxwell system we introduce a 
generating set of solutions, with its aid we obtain the Maxwell system in the 
form of an equation for pseudoanalytic functions of the second kind. This allows 
us to find a relation of the Maxwell system in the one-dimensional case with 
the hyperbolic pseudoanalytic function theory developed in [T^]. As a direct 
consequence of this relation we obtain an infinite system of solutions of the 
Maxwell equations. 

In the case of force-free magnetic fields we show that in fact one exact so- 
lution is sufficient to obtain a corresponding generating quartet and to write 
down the second-kind equation. This interesting phenomenon leads to an ob- 
servation regarding the quotients of solutions. Namely, we obtain a differential 
equation satisfied by the quotients of solutions. This result seems to be new 
even in the case of monogenic functions (a special case when the proportional- 
ity factor in the considered system vanishes identically). Finally, in section 5 
we show that the concept of a generating set is applicable to the Dirac system 
with electromagnetic and scalar potentials, and in this case as well the system 
can be written in another form corresponding to pseudoanalytic functions of the 
second kind. 

2 Preliminaries 
2.1 Biquaternions 

We will denote the algebra of biquaternions or complex quaternions by M{C) 
with the standard basic quaternionic units denoted by cq = 1, ei, 62 and 63. The 
complex imaginary unit is denoted by i as usual. The set of purely vectorial 
quaternions 5 = q is identified with the set of three-dimensional vectors. 

The quaternionic conjugation of a biquaternion q — Qo + (I will be denoted 

as 

q^qa-q, 

and by q* we denote the complex conjugation of q, 

q* = Rc q — i Im q. 

Sometimes the following notation for the operator of multiplication from the 
right-hand side will be used 

MPq = q-p. 

The main quaternionic differential operator introduced by Hamilton himself 
and sometimes called the Moisil-Theodoresco operator is defined on continuously 
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differentiable biquaternion- valued functions of the real variables xi, X2 and X3 
according to the rule 

3 

where dk = -fp—- 

2.2 Pseudoanalytic functions 

In this subsection we introduce some basic concepts from Bers' theory of pseu- 
doanalytic functions [2] and give a slightly varied interpretation of the notion 
of a generating pair. Precisely this different interpretation allows us to intro- 
duce the generating sets for the considered in the subsequent sections systems 
of mathematical physics. 

According to [5] a pair of arbitrary continuously differentiable with respect 
to the real variables x and y complex-valued functions F and G satisfying the 
inequality 

Im(FG) > (1) 

in a domain 17 C C is called a generating pair. This inequality means that F 
and G are independent in the sense that any complex function W defined in Q 
can be expressed in the form 

W ^ipF + iPG 

where if and ip are real- valued functions. 

For a fixed generating pair and in the case when Lp and ip are continuously 
differentiable one can define the {F, G')-derivative of the function W in the fol- 
lowing way 

W = F d^ip + G d^ij, 
where dz — -^{dx — idy). The derivative W exists iff the equality 

F dz(p + G dsiP = (2) 

holds. Here dz = ^{dx + idy). 
Introducing the notation 

FdzG - GdzF ^ _ FdzG - GdzF 
FG-FG ' FG-FG ' 

equation ^ can be written in the form of a Vekua equation 

dzW - aW ~ bW ^ 0. (3) 

Functions a and b are known as characteristic coefficients of the generating 
pair {F, G) and solutions of ([3]) are known as pseudoanalytic functions, or more 
exactly {F, G')-pseudoanalytic functions of the first kind. Solutions of the cor- 
responding equation (j2| regarded as complex-valued functions w = + iip are 
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called {F, G)-pseudoanalytic functions of the second kind. Note that by con- 
struction both generating functions F and G are {F, G)-pseudoanalytic of the 
first kind. 

Now let us consider the operator from Vekua equation — a — bC 

where by C we denote the operator of complex conjugation. Take an arbitrary 
real-valued function ip and consider the equality 

{d^~a-bC){^f)^fch^ (4) 

where / is some complex function. It is easy to see that this equality holds for 
any real- valued (/? iff / is a particular solution of ([3]). In order to be able to 
consider a general solution of Q another particular solution, say, g is needed. 
In this case we can look for solutions of (jSj in the form W = (pf + ^'9 where 
ip and ■0 are real-valued functions and / and g are particular solutions of (|3| 
if only / and g are independent in the sense explained above. Thus we arrive 
at the concept of a generating pair for a Vekua equation via the intertwining 
relation Q. We have then 

{ck-a- hC) {ipf + ipg) = fd^ip -I- gd^i^ 

which in particular gives us a relation between equations ([3]) and 



3 Generating sets of solutions for the Maxwell 
system in inhomogeneous media 

Let us consider the Maxwell equations for inhomogeneous media 

rotH = £atE+j, (5) 

rot E = -M^tH, (6) 

div(£E) = p, (7) 

div (^H) = 0. (8) 

Here e and /i are real- valued functions of coordinates, E and H are real- valued 
vector fields depending on t and spatial variables, the real-valued scalar function 
p and the real vector function j characterize the distribution of sources of the 
electromagnetic field. 

The wave propagation velocity will be denoted by c = — the refraction 

index by n = y/sji and the intrinsic impedance of the medium by Z = ■ As 
was shown in [7j, [H], introducing the notations 



grad \/c. „ grad , , , ^„ ^„ 

c=- — Z=2 — and V = ViE -I- iV^iH 

one can rewrite system ([51)- ([S]) in the form of a single biquaternionic equation 

{-dt + iD)Y - M'^Y - AP^Y* = -(^j + ^) 
c ye 
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which in a sourceless situation becomes 

{^dt + iD)V - AP'^V - AP^V* ^ 0. (9) 
Let ip he a, real-valued function. Then it is easy to see that the equality 

{-dt + iD- AP" - AP^C)[ipY] = {-dt + iD)[Lp\ ■ V (10) 
c c 

holds iff V is a solution of ©. 

Assume that {Vi,...,V6} are solutions of ([9]) independent in the sense 
that for any complex vector function V there exist real valued functions 
k = 1,2,..., 6 such that ^ ~Yl\=i'^k^k- This can be easily written as a 
condition on a corresponding determinant of a matrix formed by components of 
Vfc. Then due to (flU)) we have that V — Yl\=i V'/cVfe is a solution of © if and 
only if 

V(-at+^Z?)[¥>fc] -Vfc = 0. (11) 
^-^ c 

k = l 

This equation for real-valued functions i^fc, k = 1, 2, . . . , 6 is a Bers' equation 
for Maxwell pseudoanalytic functions of the second kind. 

Notice that in a frequently encountered in practice case fi — Const we have 
that c = Z and (j9]) turns into the equation 

{^dt + iD)V -i{V + V*)c^O, (12) 

for which it is easy to propose a triplet of independent solutions in the form 

V4 = ?ei, V5 = ie2, Ve = ies 

corresponding to a constant magnetic field. Thus, in order to rewrite Maxwell's 
system in the form ()lip it is sufficient to find another triplet of solutions. In 
some cases this can be done relatively easy. Let us consider a stratified medium, 
that is e = e{x) and hence c = ci(x)ei = 2c(x) ^i' Then the vectors 

Vi = cei, V9 = — , and V3 ~ — 
c c 

are solutions of (fT^ . Consequently, the Maxwell system in this case is equivalent 
to equation ((TT|) (for Maxwell pseudoanalytic functions of the second kind): 

1 '^1 e ^ 1 

{-dt + iD)(pi ■ cei + y^(-dt + iD)ipk ■ h y^(-dt + iD){pk ■ ieu^z = 0. 

c ^ — ' c c ^ — ' c 

k=2 k=4 

For a better understanding of this equation as well as of equation (IT^ let us 
consider solutions depending on t and x only: 

i^dt + teA)Y{t, x) - i(V{t, x) + V*(t, x))ci{x)ei = 0. 
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One can observe that equations for Vi and those for V2, V3 are not coupled. We 
have 

i^dt + ieid,)Vi - iiVi + V,*)ciei = (13) 

and 

i^dt + ieA)iV2e2 + ¥363) - ii{V2 + V;)e2 + (^3 + V^3*)e3)ciei = 0. (14) 

The first of these equations can be easily solved. Note that dtVi = and 
d,j:Vi-2ci Re Vi = 0. That is Im Vi = Const and Re T^i - ^ Re Vi = which 
gives us a general form of the component Vi in the case under consideration, 
Vi = aic{x) + ia2 where oi and 02 are arbitrary real constants. 

Now let us consider equation . It can be written as the following bicom- 
plex equation 

{-dt + ieA)^ + iciei($ + $*) = 
c 

for the bicomplex function <i> = V2 + ¥361. 

Denote by TV an antiderivative of the refraction index n and consider the 
following change of the variable x ^ = N{x). Then the function \l/(t, ^(a;)) = 
^{t,x) as a function of the variables t and ^ satisfies the following equation 

{dt + te,d^)^{t, + + ^*it, 0) = 

where C(^(a;)) = c{x). Note that introducing a new unity j = iei which ob- 
viously satisfies the equality = 1 and considering the bicomplex function in 
the form = ^I^i + \l/2j where \E'i^2 = ^1,2 + 1^1.261 with ui^2 and vi^2 being real 
valued functions we can rewrite the last equation as follows 

where ^* = 5*1 — ^2i- Finally, applying the conjugation operator and intro- 
ducing the new bicomplex function W — ^/C"^* we arrive at the equation 

where / = VC, which can be written in the form of a hypebolic Vekua equation 

d-W - jW = (15) 

where = -^{d^ — jdt) and instead of an asterisk we used bar for denoting 
the same conjugation with respect to j. This equation in the case when W has 
values in the algebra of hyperbolic numbers which we denote as Ji, that is when 
W = Wi + W2j with Wi and W2 being real valued, was introduced and studied 
in [12] and [11] in relation to the Klein-Gordon equation and in [F in relation to 
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the Zakharov-Shabat system. Notice that due to the fact that / is real valued, 
equation (jTS]) in fact consists of two separate equations for H-valued functions, 
that is it reduces to a pair of equations which were considered in the previous 
publications [S], [TT] and [T^]. In order to observe this one needs to write the 
function W in the form W = wi + W2ei where wi and W2 are 'H-valued. Then 
equation ([T5|) is equivalent to the following pair of separate equations 

dzWi - y ^ = and dzW2 ~ y ^ ^ 0. (16) 
Now gathering all the introduced transformations we have that 

and _ 

where tilde means that the corresponding original function was written as a 
function of ^, e.g., e{£,{x)) = e{x). 

Thus, all theory developed in [T^ and [TT] is applicable in this case to the 
hyperbolic pseudoanalytic functions wi and W2- As an interesting application 
let us mention a possibility to construct an infinite system of exact solutions 
of (fT5|) using the results from [12] and [11] combined with the elegant formulas 
obtained by Bers and Gelbart in the elliptic case [21. Namely, we notice that 
the pair of functions 

(F,G) = (/,j7/) (17) 

is a generating pair for both equations p6D and moreover it has a form con- 
venient for applying the results of Bers and Gelbart. Following :2! (see [Til 
Sect. 4.2] for the slightly corrected formulas), we give explicit formulas for the 
formal powers corresponding to the generating pair ([T7]) with / depending on 
one Cartesian variable ^. For simplicity we assume that zq — and F{0) — 1. 
In this case the formal powers are constructed in an elegant manner as follows. 
First, denote 

xW(e) = x(o)(e) = 1 

and for n = 1,2, ...denote 



X(")(0 = < 



n 



X(''-^\x)j4f^ for an odd n 
n / X^"^^^^ (x) {x)dx for an even n 
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^(n '^){^x)p{x)dx for an odd n 



X(^-^\x)j4^ for 



an even n 



Then for a = a' + ja" , a' , a" G M and z = ^ + tj we have 

(a, 0, z) - /(O Re (a, + Im (a, 0, z) 



where 



and 



m=0 V"^/ 

+ ( ) for an odd n 



m=0 



m— 



>{a,Q,z) - a 2^ ( '] t 



' ^"""H m ) i"*" foi' an even n. 



(19) 



For any a d H and n g N the formal power Z^"\a,0, z) is a solutfon of p^. 
Thus, the system of constructed formal powers gives an infinite system of so- 
lutions of the Maxwell equations in the case under consideration. Up to now 
it is an open question what part of the kernel of ()16p and consequently of the 
Maxwell system (fT4|) can be approximated by the obtained exact solutions. 



4 Generating solution for force-free magnetic fields 

The system describing force-free magnetic fields has the form 



rot B aB = 0, 
div B = 0, 



(20) 



where B is a vector describing the magnetic field and the commonly known 
as proportionality factor a is a real- valued function of spatial variables. This 
system is of great importance in such fields as high temperature superconductors 
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and dynamics of magnetofluids (see e.g. [1], [T3], [H]). It has been analized 
by different methods including those of quaternionic analysis (see [3], [T3] for 
the theory in the case when a is a constant and [8], [9j and [17] for some 
developments in a nonconstant case). 

Obviously system ()20p can be written in the form of a quaternionic equation 

{D + a)B = 

which is completely equivalent to (^0]) . In what follows we do not restrict 
ourselves by purely vectorial null solutions of the operator D + a and consider 
the equation 

iD + a)B = (21) 

on the class of continuously differentiable H(C)-valued functions B with the 
proportionality factor a being a complex-valued function. Applying the same 
idea as before we find out that for any scalar (complex-valued) function ip £ 
C^(n) where fl is some domain in the equality 

{D + a) {(fib) = (Dip) b (22) 

holds with b being an ]HI(C)-valued function iff 5 is a solution of (pTjl . Thus, 
chosing four independent solutions of bk, k — 0,3 we can look for a general 
solution of (PT|) in the form B — J2k=o 'fikbk where the scalar functions ip^ 
should satisfy the equation 

3 

Difik ■bk = 

k=0 

which is an analogue of the Bers equation for pseudoanalytic functions of the 
second kind. 

Nevertheless one can notice that in the case of equation (^1)) in fact it is 
sufficient to have only one particular solution b. As 6A with A being a constant 
biquaternion is also a solution of (PTjl the generating quartet bk, fc = 0, 3 can be 
proposed in the form bk = bck if only b is invertible in any point of the domain 
of interest 

Thus, we obtain the following statement. 

Proposition 1 Let b be a particular solution of i21\) . invertible in any point of 
J7 C M^. Then the general solution of 1121]) can be represented as the product 
B = b^ where the components of the M{<C) -valued function $ satisfy the equation 

3 

D<i>k ■ bck = 0. (23) 

fe=0 

Proof. Under the conditions of the proposition let us consider the function 
B = b^ — X]fc=o^*^fc^fe- Now application of the operator D + a to B gives 
{D + a)B = J2k=o ■ bck + J2k=o ■ {{D -I- a)b) Ck from where it is seen 
that B is a solution of ^ iff ^ is valid. ■ 
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Remark 2 It should be noticed that the last proposition is also valid when a is 
a full biquaternionic function. 

From the above proposition the following interesting fact about the quotients 
of solutions of (PT|) follows. 

Proposition 3 Let f and g be solutions of \21]) and f be invertible in the 
domain of interest. Then the function <i> — f^^g is a solution of the equation 

3 

D<^k • /efc = (24) 

k=0 

and if g is also invertible the inverse ^' = = g~^ f is a solution of the 

equation 

3 

^ D-^k ■ gek = 0. (25) 

fe=0 

Proof. According to the previous proposition if g is a solution of (PT|) it can 
be represented in the form g = /$ where the components of the biquaternionic 
function $ = f ^^g are solutions of Now noticing that if g is invertible as 
well the function $ is also invertible and hence consideration of / = leads 
to equation ()25p. ■ 

Remark 4 This proposition is of course also valid in the special case a = 
corresponding to quaternionic monogenic or hyperholomorphic functions. It is 
well known that in general a quotient of two monogenic functions must not be 
monogenic. Proposition\^ gives us a precise equation satisfied by the quotient. 



5 The Dirac equation 

We consider here the Dirac equation for one spin 1/2 particle under the influence 
of an electromagnetic potential, but in fact the procedure is applicable to other 
kinds of physical potentials (scalar, pseudoscalar, etc.). The Dirac equation has 
the form 

I lo9t - Y Ikdk + + 1070 + i ^ ^fc7fc ) * = 0, (26) 

\ k=l k=l / 

where 70,71,72 and 73 are the Dirac 7-matrices (see, e.g., [E]), m is the mass 
of the particle, (j) is the electric potential and Ai,A2 and ^3 are components 
of the magnetic potential The wave function is a C'-vector function 
^ = (*o, *i, *2, *3)- Equation (gHl) is considered in some domain G C M"*. 

As was shown in (see also ^ and [S]) the Dirac equation ((^ can be 
written in the following quaternionic form 

{D - dtAn + a - Ar<#"=i+™«2)$ ^ 
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where the purely vectorial quaternion a is obtained from the magnetic potential 
and the ]HI(C)-valued function $ is related to ^ by an invertible matrix 

transformation (see [13], [8] and [H]). It is worth mentioning that this form of 

the Dirac equation was recently rediscovered in [TH]. 

Consideration of solutions of (pS)) with fixed energy \E'(t,x) — ^'^(x)e*"* 

leads to the biquaternionic equation 

(L> + a + M^)W = (27) 

where W is an H(C)-valued function of three spatial variables and b = ~ + 
uj)ei — me2. 

Proposition 5 Let the biquaternionic functions Fq,Fi,F2 and be indepen- 
dent solutions of (flTp . Then W = J2t=o'fikFk where Lpo,Lpi,tp2 and (^3 are 
scalar functions is a solution of \2'T^ if and only if the following equation is 
satisfied 

3 

Y,{D^k)Fk^O. (28) 

Obviously, equation (^5]) is an analogue of equation ^ describing pseudo- 
analytic functions of the second kind. 



6 Conclusions 

In the present paper it was shown that the concept of a generating pair is not 
limited to the classical pseudoanalytic function theory and can be introduced in 
relation with a variety of first-order systems of mathematical physics. Here we 
considered the Maxwell system for inhomogeneous media, the system describ- 
ing force-free magnetic fields and the Dirac system from relativistic quantum 
mechanics. Nevertheless the approach presented in this paper and based on 
the consideration of a generating solution as an intertwining operator is clearly 
more general. The knowledge of a generating set of solutions makes it possible 
to rewrite the original system in a form analogous to the Vekua equation for 
pseudoanalytic functions of the second kind which opens the way to construc- 
tion of classes of solutions as was shown in section 3 for the Maxwell system and 
perhaps more important to the development of Bers' theory for the correspond- 
ing system of mathematical physics including such concepts as the derivative, 
antiderivative, formal powers, Taylor and Laurent series, etc. 
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